Charged perfect fluid disks as sources of Taub-NUT-type spacetimes 
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The interpretation of a family of electrovacuum stationary Taub-NUT-type fields in terms of 
finite charged perfect fluid disks is presented. The interpretation is mades by means of an "inverse 
problem" approach used to obtain disk sources of known solutions of the Einstein or Einstein- 
Maxwell equations. The diagonalization of the energy-momentum tensor of the disks is facilitated 
in this case by the fact that it can be written as an upper right triangular matrix. We find that 
the inclusion of electromagnetic fields changes significatively the different material properties of the 
disks and so we can obtain, for some values of the parameters, finite charged perfect fluid disks that 
are in agreement with all the energy conditions. 
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I. INTRODUCTION 

Several methods are known to exactly solve the Ein- 
stein and Einstein-Maxwell field equations, or to gener- 
ate new exact solutions from simple known solutions 
However, the above mentioned methods in general lead 
to solutions without a clear physical interpretation or 
to solutions that depend of many parameters without a 
clear physical meaning. Accordingly, it is of importance 
to have some appropriate procedures to obtain physical 
interpretations of these exact solutions. So, in the past 
years such procedures have been developed for static and 
stationary axially symmetric solutions in terms of thin 
and, more recently, thick disk models. 

Stationary or static axially symmetric exact solutions 
of Einstein equations describing relativistic thin disks are 
of great astrophysical importance since they can be used 
as models of certain stars, galaxies and accretion disks. 
These were first studied by Bonnor and Sackfield [2] , ob- 
taining pressureless static disks, and then by Morgan and 
Morgan, obtaining static disks with and without radial 
pressure [3, 4] . In connection with gravitational collapse, 
disks were first studied by Chamorro, Gregory, and Stew- 
art Q ■ Disks with radial tension have been also studied 
. Recently, more realistic models of thin disks and thin 
disks with halos made of perfect fluids were considered 
in [7,]. 

Several classes of exact solutions of the Einstein field 
equations corresponding to static and stationary thin 
disks have been obtained by different authors d, |^, [l^, 
M, E [S Q H H [13, with or without radial 
pressure. Except for the pressureless disks, generally the 
other disks have as source matter with azimuthal pres- 
sure (tension) different from the radial pressure (tension). 
However, in some cases these disks can be interpreted as 
the superposition of two counterrotating perfect fluids. A 
detailed study of the counterrotating model for the case 
of static thin disks is presented in [19|. 
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Disk sources for stationary axially symmetric space- 
times with magnetic fields are also of astrophysical im- 
portance mainly in the study of neutron stars, white 
dwarfs and galaxy formation. Although disks with elec- 
tric fields do not have clear astrophysical importance, 
their study may be of interest in the context of exact 
solutions. Thin disks have been discussed as sources 
for Kerr-Newman fields f20| . magnetostatic axisymmet- 
ric fields [2l[, conformastationary metrics [2^, and elec- 
trovacuum static counterrotating disks 

In all the above cases, the disks are obtained by an "in- 
verse problem" approach, called by Synge the " g-method^ 
[23 |. The method works as follows: a solution of the vac- 
uum Einstein equations is taken, such that there is a 
discontinuity in the derivatives of the metric tensor on 
the plane of the disk, and the energy-momentum tensor 
is obtained from the Einstein equations. The physical 
properties of the matter distribution are then studied by 
an analysis of the surface energy- momentum tensor so ob- 
tained. Very recently the above procedure has been gen- 
eralized in order to obtain thick disk models [1^ and also 
non-axisymmetric planar distributions of charged dust 
(26| . As we can see, by this "inverse problem" approach 
we can do physical interpretations of many known solu- 
tions of the Einstein and Einstein-Maxwell vacuum equa- 
tions, in the sense that the thin disks can act as exact 
sources for the space-time metrics given by the vacuum 
solutions. 

The purpose of this paper is the interpretation of 
a family of electrovacuum stationary Taub-NUT-type 
spacetimes in terms of finite charged perfect fluid disks. 
The paper is structured as follows. We first present, 
in Sec. II, the family of stationary Taub-NUT-type so- 
lutions of the Einstein-Maxwell equations that we will 
consider. Then, in Sec. Ill, we present a summary of 
the procedure to construct models of disks with nonzero 
radial pressure. We obtain expressions for the surface 
energy-momentum tensor and the surface current density 
of the disks. In particular we find explicit expressions for 
the surface energy density, the pressure and the surface 
charge electric density. Finally, the disks are interpreted 
in terms of a charged perfect fiuid. In Sec. IV we analyze 
the physical properties of the disks and we show that, for 
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some values of the parameters, the disks are in agreement 
with ah the energy conditions. We also study the motion 
of the disks and their stability against radial perturba- 
tions is analyzed. Finally, in Sec. V, we summarize our 
main results. 



functions and the electromagnetic potentials are obtained 
from the relations 



/ = Re £ + (8a) 



II. EINSTEIN-MAXWELL EQUATIONS AND 
TAUB-NUT-TYPE FIELDS 



The line element of a stationary axisymmetric field can 
be written in quasi-cylindrical coordinates (t,(p,r, z) in 
the form 



2*rD2j,„2 , „2A 



{dr''+dz% 
(1) 

where R, 4", W, and A are functions of r and z only. 
The vacuum Einstein-Maxwell equations, in geometrized 
units in which SttG — c = ~ Eq = are given by 



Gab — Tab 



7>ab 



= 0, 



with 



Tab — FacFu 



1 



rgabFcdF^ 



Fab — ^6, a ^ ^a.l 



(2a) 
(2b) 

(3a) 
(3b) 



For the metric ([T]), the Einstein-Maxwell equations in 
vacuum imply that R satisfies the Laplace's equation 



R, 



R^zz = 0, 



(4) 



so that the function R can be considered as the real part 
of an analytical function F{v) = R{r, z) + iZ(r^ 2), where 
1/ = r -\- iz. Thus the function F(v) defines a conformal 
transformation 



r R{r,z), 
z Z{r,z), 



(5) 



in such a way that the metric ([T]) takes the usual Weyl- 
Lewis-Papapetrou form 



ds' 



e'^{dt+Wd^Y + e-'''[R^dip^+e^'^{dR'+dZ 



(6) 

where R and Z are the Weyl's canonical coordinates. In 



this coordinates, the field equations 
the usual complex Ernst equations [27 1 



are equivalent to 
fA£ = (V5 + 2$*V$) • Vf , (7a) 



/A$ = (V£ + 2$*V$)- V$, 



(7b) 



where A and V are the standard differential operators 
in cylindrical coordinates and / = e^*. The metric 



fW,- = i?[i(Im£).c + $*4',c-***c], 
At = \/2 Re $, 



(8d) 



= V2[ii?(Im$),^+WAt,<:], 



where ^/2C, ~ R + iZ, so that V2d,(; = d,ii — id^z- 

We will consider a simple stationary Taub-NUT-type 
solution to the above system of equations given by 



* = 

A = 

W = 

At = 



I- 



-1 



c2 + 2ax + 1 



- 1 



X — y 



2kby, 
y/2ckby 



1 + c + ax 



c2 + 2ax -I- H- c2 



1 



2ax +1-1-0^ 



(9a) 

(9b) 
(9c) 

(9d) 

(9e) 



where a? + b"^ — 1 + , being c the parameter that con- 
trols the electromagmetic field. The prolate spheroidal 
coordinates, x and y, are related to the Weyl coordinates 

by 



Ri'' = k'^{x^ -l){l-y^), Z = kxy, (10) 



where 1 < a; < 00 and < y < I. This solution can 
be generated, in these coordinates, using the well-known 
complex potential formalism proposed by Ernst j27! | from 
the Taub-NUT vacuum solution Note that when 

c — this solution reduces to the Taub-NUT vacuum 
solution. 

Once the above solution is known, we can obtain a so- 
lution of the field equations Q in the original coordinates 
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by setting 



R{r, z) 


= Re F{v), 


(lla) 


*(r, z) 




(lib) 


A(r, z) 


= A{R,Z) + ln|F'(i/)|, 


(He) 


>V(r, z) 




(lid) 






(He) 


A^{r, z) 




(llf) 



where F' = dF/dv. Note also that the solution to the 
original Einstein-Maxwell equations ^ is not completely 
determined as we do not take any especific choice for 
F{v). We will do this in the next section. 



For the metric ([T]), the nonzero components of are 
5o" = [2i?2(A„- 2*,,)+ 2i?i?„ - e4*WW,, ](;3a) 

Sh - (13c) 
51 = £l^[2i?2A,,+ e4*WW,.], (13d) 

Si = ^[2i?i?,,], (13e) 

and the nonzero components of the surface current den- 
sity are 

k = -2e*-^At,,, (14a) 



III. A FAMILY OF FINITE CHARGED 
PERFECT FLUID DISKS 

In order to obtain a solution of ([2]) representing a thin 
disk at z = 0, we assume that the components of the 
metric tensor and the electromagnetic potentials are con- 
tinuous across the disk, but with first derivatives discon- 
tinuous on the plane z = 0, with discontinuity functions 

^ab 5afc,z|^^(j+ 5ab,2:|^^(,- 9ab,z\ j 

Thus, the Einstein-Maxwell equations yield an energy- 
momentum tensor T^b = T°^'° + where T^°^* = 
Qab '^(•z), and a current density Jq = jaS{z) = 
— 2e^'*"'^-'A(i,z(5(z), where S{z) is the usual Dirac func- 
tion with support on the disk. T^^™ is the electromag- 
netic tensor defined in Eq. ([5a|) , ja is the current density 
on the plane z = 0, and 

Qb = \{h'''-5l-¥-'-6t + g'^W,-g^-^X + K{f''K-r-&i)} 

is the distributional energy-momentum tensor. The 
"true" surface energy-momentum tensor (SEMT) of the 
disk, Sab^ and the "true" surface current density, j^, can 
be obtained through the relations 

Sab = jr^b^'dsn = e^-^ Qab. (12a) 

L ^ J J- "^'n = e^-* ja , (12b) 

where ds„ = y^gzz dz is the "physical measure" of length 
in the direction normal to the disk. 



= -2e*-^A^,., (14b) 

where all the quantities are evaluated at z — 0^. 

We now consider a solution with the above mentioned 
discontinuity properties by taking the Taub-NUT-type 
solution ([9]) with F{h') given by 

F{iy) = u + ax/i''^ - 1 , (15) 

where a > 0. This choice for F^v) was firstly presented 
in reference @ and leads to thin disks with nonzero radial 
pressure and of finite radius, located at z = 0, < r < 1. 
In order to obtain disks with non-unit radius, we only 
need to make the transformation r ar, where a is the 
radius of the disk. 

With this choice for F{i') the image of the disk by the 
conformal mapping ([5]) is the surface 

a^i?^ + ^a^ (16) 

so that the disks are mapped into spheroidal thin shells 
of matter and its exterior is mapped into the exterior 
of shells. We have three possible values for a: a ~ 1, 
a spherical shell, a > 1, a prolate spheroidal shell, and 
< a < 1, an oblate spheroidal shell. Therefore, the co- 
ordinates naturally adapted to the symmetry of the shells 
are, respectively, spherical, prolate spheroidal and oblate 
spheroidal. By considering the above shells as sources, 
we then seek for exterior solutions of the field equations 
d?]) and, using (fTTjl . we obtain the corresponding disk so- 
lutions in the original coordinates. Thus, for example, 
for the Taub-NUT-type metric (O, written in prolate 
spheroidal coordinates, we must choose a > 1 and the 
shell would be located ata; = a/A:>l and y = \J\ — . 
with k = \/ oP- — 1. 

Due to the fact that the metric function yV'(r, z) de- 
pends only of the prolate spheroidal coordinate y, and as 



4 



a consequence of the behavior of the derivatives of F{v), 
we can easily see that Wjz = at the disk and so we have 
= and 5? = (S^-SDW. That is, the SEMT can be 
written as an upper right triangular matrix and so can 
be easily diagonalized in terms of an orthonormal tetrad 
ea^ = {V^^M^^X^y''}, where 



ya 


= e-* (1,0,0,0), 


(17a) 








w 


= ^ (-W,1,0,0), 


(17b) 




= e*-^(0, 0,1,0), 


(17c) 




- e*-'^(0, 0,0,1). 


(17d) 



Since the vectors ea'' are the eigenvectors of the SEMT, 
the timelike vector can be interpreted as the veloc- 
ity vector of the disks and so the orthonormal tetrad is 
comoving. This is the orthonormal tetrad used by static 
observers who are at rest with respect to infinity, or "Lo- 
cally Static Observers" (LSO) f23|. 

We obtain for the metric, the SEMT and the surface 
current density the expressions 

9ab = -VaVb + WaWb + XaXb + YaYb, (18a) 



where a = a/k > 1. From (|22aP follows that the fluid has 
a barotropic equation of state. The effective Newtonian 
density, defined as g = e + 2p, is 



a(a2 + 1) + {2 + c^)a 



(23) 



The surface current density in the coordinates frame is 
given by 



Jt 



ad? + (1 + c2)(2q; + a) 



8y2fc2(l +c2)l/2 



(|)^ (24a) 



2A5(l-r^)V2,- 



(24b) 



whereas the surface electric charge density measured by 
the comoving observer is 



ac? -f (1 -t- c2)(2a a) 



4\/2fc(l + c2)i/2 



(25) 



and the azimuthal current density measured by the co- 
moving observer is equal to zero. Thus we have a family 
of finite charged perfect fluid disks with constant surface 
energy density and pressure given by (|22p . and constant 
surface electric charge density given by ([25|) . 



Sah = eVaVb + p^WaWb + PrXaXb 



ja = <jVa+jWa, 



where 



•^0 5 Pf ^ ^li Pr 



a2 



(18b) 



(18c) 



(19) 



are, respectively, the surface energy density, the az- 
imuthal pressure, and the radial pressure measured by 
this observer, and 



(20) 



are the electric charge density and the azimutal current 
density of the disk measured by this observer. 

From the expressions ([9]) for the Taub-NUT-type solu- 
tion, is easy to see that p,p = Pr = p and so the surface 
energy-momentum tensor can be cast in the perfect fluid 
form 



Sab = {e + p)VaVb + phab 



(21) 



where hab — gab — YaYb is the metric of the z = hiper- 
surface, 



a + a 
2fc2 



2a 



^a^ + 2adL -|- 1 + c2 ' 



(22a) 
(22b) 



IV. THE PHYSICAL PROPERTIES OF THE 
DISKS 

We will now analyze the physical properties of the 
disks. In first instance we consider that c = 0, so that 
(p- ^ \p- — \, In this case the electromagnetic potentials 
are equal to zero and so the disks are made of a neutral 
perfect fluid. Furthermore, is easy to see that with |a| < 1 
the energy density is always a negative quantity and so 
the disks obtained when c = never satisfy the weak en- 
ergy condition. So, in order to have disks in agreement 
with the weak energy condition, we need to consider so- 
lutions with c 7^ 0. That is, disks with nonzero surface 
electric charge density. 

From expression (|22ap is easy to see that the energy 
density is positive if we take 



a < —a < — 1. 



(26) 



Additionally, in order to have real positive expressions 
for the pressure p and the effective Newtonian density g, 
we must to impose the conditions 



2aa +l + c^ > 0, 



(27a) 



a{a^ + 1) + {2 + c^)a > 0. (27b) 
Thus, for example, if we take c = 3, a = —3 and a = 2, 
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we have that 



e = 3V2, 

P = 2V2, 

Q = 7V2, 

cr = 0.6\/l5, 



(28a) 
(28b) 
(28c) 
(28d) 



in agreement with the weak and strong energy conditions. 
Furthermore, p < e and so the disk also satisfies the 
dominant energy condition. 

In order to analyze the motion of the disks, we can 
see that the spatial components of the velocity vector V"" 
are zero with respect to the coordinates and so the disks 
are "locally static" . The motion of the disks can also be 
analyzed by considering an orthonormal frame different 
from the comoving tetrad. So, an orthonormal tetrad 
commonly used is the "Locally Nonrotating Frame" or 
"Zero Angular Momentum Observer" [29,, M\ ■ The tan- 
gential velocity of the disks measured by this observer is 
given by 



where uj 



gii(rj - uj) 

^LNRF / — n •> 

V5oi - 500511 
-.901/511 and n = V^/V^. We obtain 



b ( p\'^ \J\ — "^"^ 
2k 



(29) 



(30) 



As we can see from the above expression, with respect to 
LNRF, the particles of the disks move with superluminal 
velocities for r < rg, where 



''o 



6p2 



+ 4fc2(; 



(31) 



For instance, with the values used in 



we have that 
and so the disks will have sublu- 



0.96, 

minal velocities only in a narrow region near the edge. 
However, the above superluminal behavior of the veloc- 
ity is indeed due to the motion of the LNRF and not to 
the motion of the disks. So, we can see that for this fam- 
ily of solutions the LNRF is not a well behaved observer 



in the sense that the frame itself presents superluminal 
velocities. 

Another quantity related with the motion of the disk 
is the specific angular momentum of a particle rotating 
at a radius r, defined as h = gipaV^ ■ Thus we have 



bp 



(32) 



This quantity can be used to analyze the stability of the 
disks against radial perturbations. The condition of sta- 
bility, 



rf(/i^) 
dr 



> 0, 



(33) 



is an extension of Rayleigh criteria of stability of a fluid 
in rest in a gravitational field [11]. As we can see, the 
disks are not stable under radial perturbations. 



V. DISCUSSION 

The interpretation of a family of electrovacuum sta- 
tionary Taub-NUT-type fields in terms of charged perfect 
fiuid disks of finite extension was presented. The inter- 
pretation was done by means of an "inverse problem" 
approach used to obtain disk sources of known solutions 
of the Einstein or Einstein-Maxwell equations. The diag- 
onalization of the energy-momentum tensor of the disks 
was facilitated by the fact that it was written as an up- 
per right triangular matrix. We find that the inclusion 
of electromagnetic fields changes significatively the dif- 
ferent material properties of the disks and so, for some 
values of the parameters, finite charged perfect fluid disks 
were obtained that are in agreement with all the energy 
conditions. However, the disks are not stable under ra- 
dial perturbations. We also find that the disks have a 
barotropic equation of state. 



Acknowledgments 

The authors want to thank the finantial support from 
COLCIENCIAS, Colombia. 



[1] D. Kramer, H. Stephani, E. Merit, and M. McCallum, Ex- 
act Solutions of Einsteins's Field Equations (Cambridge 
University Press, Cambridge, England, 1980). 

[2] W. A. Bonnor and A. Sackfield, Commun. Math. Phys. 
8, 338 (1968). 

[3] T. Morgan and L. Morgan, Phys. Rev. 183, 1097 (1969). 
[4] L. Morgan and T. Morgan, Phys. Rev. D 2, 2756 (1970). 
[5] A. Chamorro, R. Gregory, and J. M. Stewart, Proc. R. 



Soc. London A413, 251 (1987). 
[6] G. A. Gonzalez and P. S. Letelier, Class. Quantum Grav. 

16, 479 (1999). 
[7] D. Vogt and P. S. Letelier, Phys. Rev. D. 68, 084010 

(2003). 

[8] D. Lynden-Bell and S. Pineault, Mon. Not. R. Astron. 

Soc. 185, 679 (1978). 
[9] P.S. Letelier and S. R. Oliveira, J. Math. Phys. 28, 165 



(1987). 

[10] J. P. S. Lcmos, Class. Quantum Grav. 6, 1219 (1989). [21] 
[11] J. P. S. Lcmos and P. S. Letelier, Class. Quantum Grav. [22] 
10, L75 (1993). 

[12] J. Bicak, D. Lynden-Bell, and J. Katz, Phys. Rev. D 47, [23] 
4334 (1993). 

[13] J. Bicak, D. Lynden-Bell, and C. Pichon, Mon. Not. R. [24] 

Astron. Soc. 265, 126 (1993). 
[14] J. Bicak and T. Ledvinka, Phys. Rev. Lett. 71, 1669 [25] 

(1993) . 

[15] J. P. S. Lemos and P. S. Letelier, Phys. Rev D 49, 5135 [26] 

(1994) . 

[16] J. P. S. Lemos and P. S. Letelier, Int. J. Mod. Phys. D [27] 
5, 53 (1996). [28] 

[17] C. Klein, Class. Quantum Grav. 14, 2267 (1997). 

[18] G. A. Gonzalez and P. S. Letelier, Phys. Rev. D 62, [29] 
064025 (2000). 

[19] G.A. Gonzalez and O. A. Espitia, Phys. Rev. D 68, [30] 
104028 (2003). 

[20] T. Ledvinka, J. Bicak, and M. Zofka, in Procceeding of 8th 
Marcel-Grossmann Meeting in General Relativity, edited 



by T. Piran (Workl Scientific, Singapore, 1999) 

P. S. Letelier, Phys. Rev. D 60, 104042 (1999). 

J. Katz, J. Bicak, and D. Lynden-Bell, Class. Quantum 

Grav. 16, 4023 (1999). 

G. Garcia R. and G. A. Gonzalez, Phys. Rev. D 69, 
124002 (2004). 

J. L. Synge, Relativity: The General Theory. (North- 
Holland, Amsterdam, 1966). 

G. A. Gonzalez and P. S. Letelier, Phys. Rev. D 69, 
044013 (2004). 

D. Vogt and P. S. Letelier, Class. Quantum Grav. 21, 
3369 (2004). 

F.J. Ernst, Phys. Rev. D 168, 1415 (1968). 

L.D. Landau and E.M. Lifshitz, Fluid 

Mec/ianics(Addisson- Wesley, Reading, MA, 1989). 

J.M. Bardeen, W.H Press, and S.A. Teukolsky, Astro- 

phys. J. 178, 347 (1972). 

S. Chandrasekar, The Mathematical Theory of Black 
Holes. (Oxford University Press, 1992). 



